tU4-7M 


MOD-ANOVAt  AN  ANALYSIS  OF 
VARIANCE  PROGRAM  FOR 
2*  AND  3*  FACTORIAL  EXnBRIMENTS 


DISCLAIMER 


The  findings  in  this  document  are  not  to  be  construed 
as  an  official  Department  of  the  Army  position  unless  so 
designated  by  other  authorized  documents.  Comments  or 
suggestions  should  be  addressed  to: 


Commander 

USA  Concepts  Analysis  Agency 
ATTN:  Director  of  Methodology,  Resources  and 

Computation 
8120  Woodmont  Avenue 
Bethesda,  Maryland  20014 


UNCLASSIFIED 


irCU^lTY  CL  A^SlFlCATfON  OF  TM'5  •AOt  'Wh»n  fmtm  Fnfm’t^d) 


REPORT  DOCUMENTATION  PAGE 


RE,F*0£T  NUMPF  « 


3App  ' RhAD  fNSTRL’CTIONS 

COMPLETING  FORM 

S GOVT  ACCF.SSION  NO^  « fc  CJP I F- S t • s C A T A ^OG  N U M y E R 


' CAA-3-77-2  _ ; 

4.  TiTLfc  

An  Analysis  of  Variance  Program 
2^ -Hd  ^^actorial  Fxperiment<^  y 


au.rMoiy.i 

Carl  B./^Bates  '' 
Jerry/  Thomas 


! Do7:umentation 

fc  PrRFOBMING  ORG  REPORT  NUMBER 
I • "cVl^R  ACT’oR  grant  NUMBERfaj 


9 PERFORMING  ORGANiZATtON  NAME  AnO  ADDRESS 

US  Army  Concepts  Analysis  Agency 
8120  Woodmont  Avenue 
Bethesda , Maryland  20014 

n.  CONTROLLING  OFFICE  NAME  AnO  ADDRESS 


<0  PROGRAM  element.  PROJECT.  TASK 
AREA  ft  WORK  UNIT  NUMBERS 


7997B 


I iz  report  oats 


US  Army  Concepts  Analysis  Agency  j ^ /■  W77 

8120  Woodmont  Avenue  **“**^P'^” 

Bethesda,  Maryland  20014 J-Z^' 3 - . 38 

u.  MOnTto^NC.  A<TenCY  name  ft  ADORESS/rf  dlUfttmtU  from  C.^rtrotting  Otflro)  HtT'  CLASS,  (of  fh/«  roport) 


^ -Ln.  KUMMfiOF  PAG« 

OJz^3^J:'  38 


Unclassified 


15«  DECLASSIFICATION  DOWNGRADING 
SCHEDULE 


16  DISTRIBUTION  STATEMENT  (of  thl»  Roport) 


{Approved  for  public  release;  distribution  unlimited. 


17,  DISTRIBUTION  STATEMENT  (of  th*  mbttemct  •n/er#d /n  R/r»c*  70.  If  from  Report) 


1*  supplementary  notes 


I 1^  KEY  WORDS  fConfInii#  on  r*v*r*#  tidm  If  nocmmamry  mnd  Irlontity  by  block  numbtr) 


2 Lc  r,-tK  .rvefr  tv  i Z i' 


Computer  program  Aliases 

Analysis  of  variance  Modular  arithmetic 

Fractional  factorial  designs 

Factorial  effects  o,.i  *1  „ /-i 

Confounding L L=.  r,  t K ,cveFr  tv  \ ' 

20-  AVSYRACT  fConlfftum  m\  rororM  vlWR  If  noK^mmfr  «»<#  M#nf/ff  fty  block  numbor)  - “ ' ■ 

Modulo  Analysis  of  Variance  (MOD-ANOVA)  is  a computer  program  for(^'^_and  3^  full 
or  fractional  factorial  experiments.  The  program  applies  modular  arithmetic  to 
cell  identification  numbers  to  classify  and  sum  the  observations.  MOD-ANOVA  is 
fast  and  efficient  and  is  extremely  simple  in  structure.  The  complete  program 
contains  less  than  130  lines  of  code.  The  program  is  written  in  FORTRAN  V and  is 
operational  on  the  UNIVAC  1108  computer,  but  it  should  be  transferable  to  any 
computer  compatible  with  FORTRAN. 


DO 


ftWTIONOf  • WOW  «»  IS  OBSOL  FTP  UNCLASSIFIED  ' 

^ ^ SPCUPITY  CLPSSIPICATioW  op  this  P^CE  (WSw,  P.I.  Enl»r»rf) 


DOCUMENTATION 

CAA-D-77-2 


MOD-ANOVA:  AN  ANALYSIS  OF  VARIANCE  PROGRAM  FOR  2'^  AND  3"  ' 

FACTORIAL  EXPERIMENTS 


JULY  1977 


Prepared  by 

methodology,  RESOURCES  AND  CuMPuTAl ION  DIRECTORATE 

US  ARMY  CONCEPTS  ANALYSIS  AGENCY 
8120  WOOOMONT  AVENUE 
BETHESDA,  MD  20014 


CAA-D-77-2 

MOD-ANOVA;  An  Analysis  of  Variance  Program 
for  2^  and  3^  Factorial  Experiments 
ABSTRACT 

Modulo  Analysis  of  Variance  (MOD-ANOVA)  is  a computer  program  for  the 
analysis  of  2'^  or  3^  factorial  experiments.  The  program  is  applica- 
ble for  either  full  or  fractional  factorial  designs.  The  program 
applies  modular  arithmetic  to  the  cell  identification  numbers  to 
classify  and  sum  observations.  MOD-ANOVA  is  fast,  efficient,  and 
extremely  simple  in  structure.  The  complete  program  contains  less 
than  130  lines  of  code.  Th^-  program,  written  in  FORTRAN  V,  is  opera- 
tional on  the  UNIVAC  1108  computer,  but  it  should  be  transferable  to 
any  computer  compatible  with  FORTRAN  V. 
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1.  INTRODUCTION 

Modulo  Analysis  of  Variance  (MOD-ANOVA)  is  a computer  program  for 
performing  the  analysis  of  variance  computations  for  a fractional 
factorial  design.  The  program  is  applicaole  for  2^^  and  3"  designs. 

It  is,  to  the  knowledge  of  the  authors,  unique  in  its  application  of 
modular  arithmetic.  The  program  originated  from  a requirement  for 
the  analysis  of  variance  of  a one-ninth  fractional  factorial  of  a 3^ 
design.  Although  the  program  was  developed  for  fractional  3^  de- 
signs, it  is  also  applicaole  to  fractional  2^  designs  and  full  2^  and 
3^  designs. 

A common  procedure  for  performing  ANOVA  computations  for  frac- 
tional 2^  and  3^  designs  is  to  perform  an  ANOVA  on  a lower  order  full 
design  and  then  use  the  aliases  to  associate  sum  of  squares  with  the 
proper  factorial  effects.  This  works  well  for  2^’  fractional  designs, 
but  the  procedure  becomes  extremely  tedious  for  fractional  designs 
as  n becomes  large.  An  alternative  method  for  performing  the  compu- 
tations is  to  employ  the  General  Linear  Hypothesis  Model  and  apply 
least  squares.  For  large  designs  this  involves  the  construction  of 
extensive  design  matrices  and  the  inversion  of  many  high  order  mat- 
rices. Because  of  the  above  cited  criticism  of  existing  computa- 
tional methods,  an  alternative  computational  scheme  was  sought.  This 
report  documents  the  computational  procedure  ultimately  developed. 
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Advantages  of  the  procedure  are  its  simplicity  and  computational  ef- 
ficiency. 

2.  BACKGROUND 

Theory  of  fractional  factorials  for  2^  and  3*^  designs  was  devel- 
oped by  Finney  (1945)  and  (1946).  The  interested  reader  is  referred 
to  the  more  readily  available  references  of  Kemothorne  (1952),  Davies 
(I960),  Winer  (1962),  or  Kirk  (1968).  Kempthorne  is  more  theoretical 
than  the  other  references  and  uses  finite  group  theory  in  the  devel- 
opment of  fractional  factorials.  Cochran  and  Cox  (1957),  Connor  and 
Zelen  (1957)  and  (1959),  and  Davies  (1960)  contain  fractional  facto- 
rial designs  for  2^  and  3*^  experiments.  (For  References,  See  Appendix  B.) 

3.  NOTATION  AND  DEFINITIONS 

Consider  an  m*^  factorial  experiment,  where  n is  the  number  of  fac- 
tors and  m is  the  number  of  levels  of  each  factor  and  is  restricted 
to  either  2 or  3.  Denote  the  general  cell  identification  by  lower 
case  letters  abc  ...,  and  denote  factorial  effects  by  upper  case  let- 
ters A,B,C, For  2^  designs,  lower  case  letters  may  take  on  0 or 

1;  for  3"^  designs,  they  may  take  on  0,1,  or  2.  A fractional  repli- 
cate involving  a subset  of  m^"'"  factor  level  combinations  of  the  full 
design  is  termed  a (l/m'")  replicate.  The  total  degrees  of  freedom  of 
a full  design  is  m"  - 1;  the  total  degrees  of  freedom  of  a d/m'") 
replicate  design  is  (m^  - D/m*". 

Let  I be  the  usual  identity  element.  The  defining  contrast  of  a 
(1/2'^)  replicate  consists  of  I equated  to  (2*^  - 1)  groups  of  letters 
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“11  “21  “il 

A B C where  a .j  ,2j  ,a3j  . . . take  on  the  values  of  0 or  1 

and  aU  = gU  = Pq^  g (l/B"^)  replicate,  I is  equated  to 

(3'^  - l)/2  groups  of  letters  and  ’^2j ’^3j  ‘ ‘ ^ 

Connor  and  Zelen  (1957)  and  (19bb).)  The  conventional  modular  nota- 
tion is  used  for  expressing  interaction  components.  Modular  notation 
is  used  Dy  Kempthorne  (1952),  Winer  (1962),  Hicks  (1964),  and  Kirk 
(196B)  and  readily  lends  itself  to  extension  and  generalization. 

Winer  (1962)  and  Kirk  (196B)  give  equivalences  between  modular  nota- 
tion and  the  older  Yates  (1937)  notation.  For  the  reader  familiar 
with  the  Yates  notation,  the  following  equivalences  are  given. 

Modular  notation  Yates  notation 


AB 

AB(J) 

AB'^ 

AB(I) 

ABC 

ABC(Z) 

ABC^ 

ADC(Y) 

AB^C 

ABC(X) 

AB^C^ 

ABC(W) 

Conventionally  the  exponent  of  the  first  letter  in  the  group  of 
letters  is  always  made  unity  since  it  can  be  shown,  for  example, 
A^BC  = (A^BC)^  = A'^B^C^  = AB^C^  modulo  3.  Each  component  of  the  3*^ 
designs  has  two  degrees  of  freedom.  Therefore,  (m'^  - l)/2m'^  compo- 
nents can  be  partitioned  from  the  total  sum  of  squares  of  a (1/3'') 
repl icate. 
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4.  COMPUTATIONAL  PROCEDURE 

The  computational  procedure  of  MOD-ANOVA  is  based  on  a linear 
function  comparable  to  Kempthorne's  rule  (1952)  for  confounding. 
Effects  or  components  of  effects  in  an  ANOVA  model  (of  a full  design 
or  after  fractionating)  are  identified.  The  program  then  computes 
the  sum  of  squares  for  each  effect  identified.  For  m = 2,  the  sums 
of  squares  are  for  whole  effects;  for  m = 3,  the  sums  of  squares  are 
for  whole  main  effects  or  interaction  components.  In  the  latter 
case,  each  sum  of  squares  has  two  degrees  of  freedom. 

Consider  the  linear  function 

E “ ^2j^2  ^nj^n’  " 1.2,. ..,h,  (1) 

where  the  a^j  coefficients  are  as  defined  in  Section  3 above  and  h is 
the  number  of  effects  to  be  calculated  for  the  ANOVA.  The  variables 
(x^.i  = 1,2,... ,n)  represent  the  ABC...  also  defined  in  Section  3 
above.  The  L-function  is  applied  to  each  abc....  cell  of  the  design 
for  each  effect  to  be  calculated.  This  results  in  a partitioning  of 
the  cells  into  m parts,  each  part  having  an  equal  number  of  cells. 
That  is,  evaluate 

L = b (mod  m);  b =0,1,...,  m-1.  (2) 
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For  example,  if  m = 2 we  have 

L = 0 (mod  2) 
and 

L = 1 (mod  2 ) . 

If  m = 3,  we  have 

L = 0 ( mod  3 ) , 

L = 1 (mod  3 ) , 
and 

L = 2 (mod  3). 

The  observations  within  each  of  the  m partitions  are  then 

summed,  squared,  and  divided  by  the  number  of  the  observations  which 
went  into  the  sum.  For  example,  consider  the  ABC^  component.  Then, 

L = Ixj^  + 1x2  + 2x3.  For  L = 0,  1,  and  2,  denote  the  three  sums  by 
(ABC^)q, (ABC^)^ , and  (ABC^)2,  respectively.  Denote  tne  three  sums 
squared  divided  by  the  number  of  observations  by  (ABC^)g,  (ABC^)j, 
and  (ABC^)^.  The  sum  of  squares  due  to  the  interaction  component 
ABC^  is  then 

SSlABC-^l  = (ABc2)g  f (ABC^lJ  + (ABC^)^  - CF , (3) 

where  CF  is  the  usual  correction  factor.  The  sum  of  squares  due  to 
the  ABC  interaction  is  the  sum  of  the  four  components'  sum  of  squares 

SSlABCl  = SSlABCl  + SSlABC^j  + SSlAB-^Cl  + SSIAB^C^I  (4) 


5 


CAA-D-77-2 

The  above  described  computations  are  illustrated  in  the  following 
numerical  example, 
b.  NUMERICAL  EXAMPLE 

a.  Problem  Description.  Consider  a one-third  replicate  of  a 3^ 
factorial  experiment.  The  243  cell  identifications  of  a full  design 
are  shown  in  Table  1.  Suppose  the  fourth  order  interaction  (ABCDE) 
is  selected  for  confounding.  Selecting  the  principal  olock  gives  the 
resulting  one-third  replicate  defined  by  the  81  circled  cell  identi- 
fications. (See  Connor  and  Zelen  (1959).)  Note  that  the  sum  of  the 
cell  identifications  equal  0 (mod  3).  The  data  for  the  numerical 
example  are  given  in  Table  2. 

The  aliases  are  obtained  by  using  the  defining  contrast  I = ABCDE 
and  the  rule  of  generalized  interaction  for  3^  systems.  Aliases  of 
the  main  effects  and  the  first  order  interaction  components  are  given 
in  Table  3.  The  Effect  Column  provides  the  needed  coefficients  for 
the  L-function  of  equation  (1).  Each  of  the  effects  defines  an  L- 
function.  Therefore,  we  have  25  L-functions  enumerated  in  Table  4 
below.  The  coefficients  are  the  exponents  of  the  Effects. 


■P. 
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Table  1. 


Full  Replicate  of  a 3^  Factorial  Experiment 


Levels  of  the 
factors 

o 

o 

t’l 

02  ] 

Bo 

tl 

B2 

Bo 

n 

Be 

*0 

Co 

C2 

iloooi 

OOTDO 

00200 

00001 

00101 

050^ 

00010 

s 

BSisin 

^Qg 

Cooof|) 

O0TT2 

00212 

i* 

c° 

01000 

OiiOO 

(^12^ 

OiiiOl 

i^Too^ 

0lT02 

01202 

giQio 

<ml5> 

OTTTO 

(gfoit) 

oTTTi 

01211 

01012 

0iU2 

01220 

01021 

(glMj 

01Q22 

*^^{72 

B2 

i? 

C2 

s 

o?Tfn 

02201 

02002 

02122 

(g2^ 

(2.201^ 

0?TT0 

02210 

02011 

02in 

(^2211) 

li 

02020 

<^222^ 

m 

(tS202t> 

07T72 

02222 

''l 

Sq 

C° 

C2 

10000 

10100 

(J02^ 

IQQOl 

(Jooo?) 
it: '^2 
10202 

lOQlO 

(IjoTI^) 

iTTTo 

UTTTl 

10211 

10012 

igm 

do^ 

10220 

10021 

10121 

(J022J) 

($0^ 

10772 

C2 

UQQO 

(JIlM 

iTToo 

dlo^ 

i'llTl 

11201 

11002 

11102 

CQ2^ 

IS 

11210 

11011 

mil 

(Xltti) 

fl 

11020 

1II20 

(11225) 

li? 

(0025) 

iTr72 

11222 

C° 

c.^ 

CX200> 

irnio 

12200 

12001 

'422PI> 

12010 

ILllg 

<iplj> 

12212 

122tg 

2^ 

12221 

12022 

$2^1) 

Bo 

c° 

C2 

<Jo^ 

20201 

20002 
2:  I )2 
Cq20> 

Cobv^ 

Trno 

20210 

20011 

20111 

(^02T> 

?sai2 

(?011^ 

2T7T2 

20020 

20120 

($0^ 

20222 

Bl 

C° 

c]  _ 

2TT30 

21200 

21001 
2UT1 
'1 120 1 

21202 

21010 

2U',0 

CT2'i> 

21Q11 

<mtj) 

2T7T1 

(?ToT^ 

2TTI2 

21212 

?iZ20 

($1125) 

2T7T0 

(fl025) 

2Tt!i 

21221 

21022 

niii 

C$J^ 

Bp 

C° 

r' 

C2 

22000 

22100 

f?22r'^ 

22001 

<^210]) 

2?7m 

22202 

2.210 

22710 

27TT1 

22211 

22012 

2iLU2 

($2^ 

22220 

22021 

($2^ 
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Table  2.  One-third  Replicate  of  a 3^ 


Levels  of  the 
factors 

Bq 

t . _J 

Co 

Cl 

^2 

■a 

El 

E2 

wm 

El 

^0 

Bo 

C° 

1996 

2016 

1978 

2398 

2088 

f^O 

Cl 

2636 

1978 

1818 

2186 

1806 

2838 

2382 

2499 

&2 

2460 

2C76 

L 

' 

279? 

2115 

2649 

2655 

2003 

Bo 

Co 

Cl 

2384 

2062 



1687 

2304 

1681 

2581 

2122 

2327 

Bl 

Co 

Cl 

C2 

2012 

■ 

1942 

2368 

«2 

Co 

C2 

1651 

2195 

1864 

2391 

2064 

1634 

2388 

1665 

1685 

1 

Bq 

Co 

C2 

1877 

1511 

■ 

2379 

2077 

2462 

1882 

Bl 

■ 

1562 

2296 

19R? 

2308 

2082 

1776 

2280 

2080 

1652 

B2 

Co 

Cl 

JiL^ 

2165 

2096 

1595 

2061 

1501 

2007 

I860 

1625 

1889 

I 
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Effect 

I = ABCDE 

A1 iases 

ABCDE 

A^B^C^D 

A 

Ab2c2o2e;2 

BCOE 

8 

AB^COE 

ACDE 

C 

ABC^DE 

ABDE 

0 

ABCD^E 

ABCE 

E 

ABCDE^ 

ABCO 

AB 

ABC^O^E^ 

CDE 

Ab2 

AC^O^E^ 

BC^D^E^ 

AC 

ab^cd^e^ 

BDE 

AC  2 

AB^D^E 

BC^DE 

AD 

AB^C^de^ 

BCE 

ad2 

AB^C^E^ 

BCD^E 

AE 

AB^C^D^E 

BCD 

Ae2 

AB^C^O^ 

BCOE  2 

BC 

Ab^C^DE 

ADE 

BC^ 

AB^DE 

AC^DE 

BD 

AB^CO^E 

ACE 

BD^ 

AB^CE 

ACD^E 

BE 

AB^CDE^ 

ACD 

BE2 

AB^CD 

ACDE^ 

CD 

ABC^O^E 

CD 

< 

CO^ 

ABC^E 

ABD^E 

CE 

ABC^DE^ 

ABO 

Ce2 

ABC^D 

ABDE^ 

DE 

ABCD^E^ 

ABC 

DE^ 

AHCD^ 

ABCE^ 
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Table  4.  Coefficients  of  L-Functions 


Effect 

&2 

^3 

34 

^5 

Effect 

d2 

^3 

34 

^5 

A 

1 

0 

0 

0 

0 

BC 

0 

1 

1 

0 

0 

B 

0 

1 

0 

0 

0 

BC^ 

0 

1 

2 

U 

0 

C 

0 

0 

1 

0 

0 

BD 

0 

1 

Q 

1 

0 

D 

U 

u 

0 

1 

u 

BD2 

u 

1 

U 

2 

U 

E 

0 

u 

0 

u 

1 

BE 

0 

1 

0 

0 

1 

AB 

1 

1 

u 

0 

u 

BE2 

0 

1 

0 

u 

2 

AB2 

1 

2 

u 

0 

0 

CD 

u 

u 

1 

1 

0 

AC 

1 

0 

1 

0 

0 

CD2 

0 

0 

1 

2 

0 

AC2 

i 

0 

2 

u 

0 

CE 

u 

Q 

1 

0 

1 

AD 

1 

0 

0 

1 

u 

CE2 

u 

0 

1 

0 

2 

AD2 

i 

0 

0 

2 

u 

DE 

0 

0 

0 

1 

1 

AE 

1 

u 

u 

0 

1 

DE^ 

0 

0 

0 

1 

2 

AE2 

1 

0 

0 

U 

2 
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b.  Program  Input.  The  program  input  consists  of  a read  format 
for  reading  cell  ID  and  data,  a format  for  writing  a-jj's,  b (mod  m) 
sums  and  the  sums  of  squares,  the  coefficients  of  the  L-function, 
cell  identifications,  and  the  data.  Input  format  descriptions  are 
given  below.  Following  the  description  is  an  illustration  of  the 
program  input  for  the  above  described  example. 


Card  Type  Purpose  Format 

1 Format  to  read  cell  ID  and  data  (10A6) 

2 Format  to  write  a^j's,  b (mod  m)  sums  (10A6) 

and  sums  of  squares 

3 Read  number  of  factors,  number  of  levels,  (415) 


fractionating  (m*^),  and  number  of  effects 
to  be  calculated  in  addition  to  the  main 
and  two  factor  interactions 

4  Cell  ID  and  data  values  (According  to 

format  on 
Card  Type  1) 
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Numerical  Example  Input 

@RUN,/TP  A106C.F3880A7997B, UNCLASSIFIED, 2, 200 
OASG.A  06M0DAN0VA. 

OXQT  06M0DAN0VA.VAR 
(1X,SI1,F9.0) 

(1X,5II,4X,4F14.4) 
b 3 3 

00000  1998. 

00012  1638. 

00021  1859. 

. (See  output  for  complete  data) 


22221  1625. 

•P 


I 

» 
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c.  Program  Output.  Program  output  consists  of  the  input  data, 
coefficients  of  the  L-functions,  and  their  sums  of  squares,  along 
with  the  sum  of  squares  of  all  observations  and  the  correction  fac- 
tor. The  difference  of  the  latter  two  gives  the  total  sum  of 
squares.  Output  from  the  above  numerical  example  follows. 
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atS6,*  0*H0ri»N0V*> 


0>OT  U4M00 *N0 V 4 • V t a 


OoOuO 

1996, 

Ooo  1 ^ 

1438< 

QQOi  1 

1 85’' 

OU  1 

1978 

OUl  1 1 

2088 

00  I 20 

2 197, 

0U2U  k 

20  18 

Ou?  1 u 

2398, 

0U222 

2330' 

U 1 002 

18  18 

0101  t 

1 8Q8  1 

0 |02u 

1 97  2. 

0 1 Ui  1 

l«78, 

U 1 1 1 u 

2 105, 

0 112* 

2‘<9V  , 

U 1 ?UC 

2635. 

0 12  1* 

2830. 

72  1 

2 382  . 

U2O0  1 

2070. 

0 70  l U 

2115. 

0*022 

2003  . 

02  1 UO 

2‘*60. 

02  1 1 2 

259  1. 

02  1 2 1 

2655. 

02202 

2792  . 

022  1 1 

268  V. 

02  2 20 

2S62  . 

1 UOo* 

1 68  7. 

igoi  1 

1681. 

1 O02U 

1 822  . 

1 0 1 0 1 

2062. 

1 0 1 10 

2 308. 

10  12* 

2 32  7, 

i 020u 

7308  , 

1 02  1 * 

256  1 . 

1 022  1 

2 122. 

1 lOol 

1 62  3. 

1 1 0 1 0 

1881, 

1 1 022 

1886, 

1 1 1 oO 

20  12. 

11112 

1982, 

1 1 1 2 1 

2366. 

1 1 202 

2203, 

II  2 1 1 

2385, 

1 1220 

__J  9 9 6 , 

1 2000 

1661  , 

12012 

1 638, 

12021 

1666. 

14 


I2l0i 

ia6<t 

mi) 

20»'« 

12I2U 

23ttl 

122U1 

21V5 

12210 

2391 

- ill 42 

.. 

2000  1 

ISl  1 

20010 

1*17 

20022 

iat2 

20100 

1177 

201  12 

2379 

20121 

29*2 

202U2 

2oai 

2021  1 

22ea 

20220 

2077 

210U0 

1S*2 

21012 

177* 

21021 

2080 

21102 

1 982 

21111 

2082 

21120 

2280 

21201 

229* 

2 1 2 lU 

2308 

21222 

1*62 

22002 

lS9b 

2201  1 

iSOl 

22020 

1 8*0 

22101 

209* 

22  110 

20S  1 

22122 

1889 

22200 

2 1*8 

2221  2 

2007 

22221 

1*28 

CA,''-D-77-2 


niMriD  Ml  i(mju  Ml  7iMon  mi  sum  Of  sau»HES 


Ancnt 

MF  AN 

1 jOju 

1 6 .0000 

bis  1 3 

u 1 ouu 

SS  7 1 *1 .0000 

bb<7H7 

Uu  1 00 

. JOUU 

b90bU 

U JO  1 u 

ShmBS.uOQO 

b7  1 7o 

OOOu  i 

S7U76.U000 

bbb  7b 

1 1 0 JO 

S‘*0‘*b.OOUC 

bb772 

1 20UU 

b32Ub.u0J0 

bH  8 9 1 

t u 1 00 

b62d3 . 0000 

bs  320 

1 O?u0 

Sb 1 u 3 . OOuO 

b6  79Q 

1 UO  1 0 

bbfttt  7 .0000 

bb7«6 

1 00^  J 

bfc2 JV.OOUO 

b62Q‘> 

1 uou  1 

bbil6Jt0000 

bbbs  1 

1 OOUi 

bb /HU .0000 

b72  1 7 

u II  uJ 

b6 366 . 0000 

bbSQV 

0 1 2UU 

bbl62.0Uo0 

b63B'» 

0 1 0 1 0 

bb'^  1 U .0000 

bb38b 

0 1 020 

bbU‘<7.0U00 

b6Q90 

□ 1001 

SbbtM.OOOU 

bS830 

0 1002 

bb /b2 .0000 

b6  1 00 

001  10 

SS3»y ,OOQQ 

‘»2S0'< 

00  120 

b 36  36 . 0000 

b6833 

OOlUl 

bbU6b . 0000 

5b723 

0 0 1 Ij  2 

b6b3o . 0000 

bb29  1 

0001  1 

b72H6.oOOO 

bS9fl6 

Qno  1 i 

bMRVV.Jooo 

b6SBtt 

R 

TOTAl 


Bf  1 N 


t i*»79*QQ? 


ncoo 

b 7 Jb I . ooon 

1 1 3 7 S40 • 2222 

QOOO 

b62 1 9.0000 

s7a2tQ7*ti 

OOuO 

60903. ooon 

36S8b3H. 7S07 

U JOO 

b622b • 0000 

1 3 7sOb  t bbb6 

0000 

bs.22V.U00n 

7 1 408.9630 

oooo 

blt063.UU00 

3J0932.bl8b 

no  jo 

bS7as.UuuO 

6J0s3b.6296 

OUOO 

b727  7 .0000 

1 6771?t  UbZ 

0000 

bb  V8  7 . ooon 

b27S3.629A 

0000 

b6S0  7 . 0000 

1 12a2.00UQ 

nouO 

bbS32 .ooon 

1 bbOS  >6667 

OOUO 

b6S79.00U0 

1 68  SB • 963U 

0000 

9<4  8 7 b . 0000 

1 032 1 6 . 9630 

0000 

bb7Q5.0Q00 

93S7»ASiS 

OOOU 

bfc  32  9 , OOUO 

3bSSH.6667 

OOUO 

b708b.UU0O 

70329.0711 

OOUO 

b*7S3.0U00 

bS20b.Bbl9 

uouc 

b7  Sd6 .uUOO 

S 39387.8S19 

0000 

b602« . 0000 

2S99.bbb6 

QUQU 

S7U&2.UaUQ 

7 32^16^7 .6296 

OOUO 

b7S 1 1 .UOOO 

3062S0«9430 

0000 

b629o . OOUO 

6S 1 2>94JO 

0000 

b60bV . 0000 

28972.6667 

QOUO 

bSbsa.OUOD 

99993.1  ai>2 

oooc 

b*39J.Qoon 

bis73.85l9 

J3J3 

. Jb«3*U9 
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Output  from  the  program  gives  the  required  sums  of  squares  to  con- 
struct the  desired  analysis  of  variance  table.  For  example,  sums  of 
squares  of  whole  effects  equals  the  sum  of  the  sums  of  squares  of  its 
components.  The  four-degree  of  freedom  whole  effect  of  the  AB  inter- 
action is  the  sum  of  the  two  two-degree  of  freedom  components  AB  and 
AB^.  That  is. 


SS(AB|  = SS|AB|  + SSlAB^l. 

The  corrected  sum  of  squares  of  total  is 

SS[T  I = TOTAL  - MEAN. 

The  sum  of  squares  of  the  residual  is  then 

SS(R  I = SS(T  I - Y (SS  of  all  fitted  effects). 

The  desired  analysis  of  variance  for  the  one-third  replicate  of 
the  3^  fractional  factorial  design  is  shown  in  Table  5 below. 
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Table  5.  ANUVA  for  the  (1/3)  x 3^  Design 


Source  of  Degrees  of  Sum  of  Mean 

Variation  Freedom  Squares  Square 


F -ratio 


9 


PROGRAM  FLOW  CHART 


CAA-D-77-2 


CAA-D-77-2 


CAA-D-77-2 


f Subroutine  It  > 
^Generates  L-functior 
V Coefficients  J 


Generate  L-function 
Coefficients  for  Main 
Effects 


Generate  L-function 
Coefficients  for  ''^wo 
Factor  interactions 


Count  Number  of  Sets  of 
Coefficients  Generated 
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7.  PROGRAM  LISTING 


O^MOOkMOVt-  « >B 

» ^I"Ptr  kNkLYITS  or  v«o!*NCr  PBOOPk"  IJ^ING  HOOOLO  kRITMUrTIC 

c 

nnuBir  opre  n;  J 31,  .cr.  ts.b 

1 NTf  GTB  » . » . « 


OT•<f^rIO^  .rrik^rt.n  (»7M.?nSt*7S>  .7ISI«Tr».77';(»Ti;) 

tnxfNrTON  sil'it  .TS(«7S1 

n*T»  «,/•»•. •a*. 
or  »oi  a .1 ) »i  ( T ) .1  r 1 . m 1 
or  »D  I*  . 1 I ( f p 1 7 I ,I  1 . 1 II I 


10« 

\ 

rOB^kT  1 Ju*f  1 

M« 

C*T* 

A /4  7A,|,|)  t 

\2m 

o ; 

0 

1 5. 

c 

IH* 

c 

Su»(7ro  or  r n roos  r i^r  . rokcnoN 

IS* 

c 

U* 

or  «0  (k  . 7 1 *( .«  . Nr  .Nf  4 

IM 

« 

roo»i4T  (Oia  1 

It* 

c 

M* 

c 

trNCPATF 

t-FUNCTION  C OFf*"  IC  Tf 

20* 

c' 

21* 

CALL 

7*  IN  ♦ A . M) 

??• 

N 77 

- M ♦ \ 

? *• 

H H 

♦ N€* 

2 

Nf  -- 

(M«  1 /^r 

2S* 

Ml  : 

Nt  /•< 

?fc* 

c 

2 /• 

c 

(MfCK  ^’Ofc 

anC  PPAC  AOniTT^NAL  FrrrcTS 

28* 

c 

29* 

T*- 1 

NT^A.^'Q.  O)  f»0  TO  \\7 

!0* 

DO  1 

T » I ' N 7 ’ 

n* 

or»o,c,.,,,»,j,j,,jr).m 

12* 

4 

F00M4T  ( 1 X . Ui  1 1 1 

51* 

1 

C 0N7 I NUf 

^ M • 

* !? 

no 

so  1 : i.Nf 

3S* 

c 

U* 

c 

Pt*D 

rrt  L 

rn  ANO  DATA  VAlUf 

5 ?• 

c 

1ft* 

Pf AO f ^ •CHT » « I M 7| • T 7 =l*N»,y 

19* 

•RI Tt I S ,f “T M X C I7 » . T y =l,N».y 

TSS 

• TSS  ♦ y*r 

23 


• j* 

4 ?• 

« 4« 


COMPOTE  '''''' 

ro  J - 

t ' U 


no 

t t* 


199  CON^l*^^*^ 
. - -nn IL 


4 9* 

,r  ,U-1'V0.?".’ 

in 

70. J>  - ''''■'’' 

ftl* 

GO  TO  701' 

<>?• 

2n 

7J.J1  - 71  '•" 

c»7* 

GO  TO  7uo 

P 4 • 

»n 

72. J.  - 77. J' 

SS- 

2 0^ 

CONTINUE 

Sf  • 
S 7* 

sn 

CONTINUE 

S9* 

t CO.PU’f 

tn* 

n« 

f 

fe 

fe<*« 
fcS* 
fc(!» 
f '• 
k-  ft* 
t 3* 

n* 

71* 

72* 

7 '• 
7»* 
7S« 
7t* 
77» 

7 9» 
79* 
AT* 
ftl* 

R 2* 
B '• 
ft*** 
»S* 
B6* 


7S  ^ 7 

CP  ?<.-.  - 

gPTT»  ' 


. I 7' 


) 

T7»  ' c . T . T 1 ' 

,aiTr  .(■  .'1  isin.T 

rnP"**’ '”  •'■‘‘,1', . •■>1-1''  “ 

r 0P»*  * ^ I I • ' 

* r /X  inn  j 


• ^7 


.,,Mor 


.,,-00 


0* 


T <.  (J  ) = ' ' ', 

P - . I-.X' 

roNTiNur 

c,  r OB"*  ^ 

1 


<c  r » 


5B«*  T I / I*  • 7^  f n . T<  > 

on  lOJ  J n .''  '■*” 

J ,gP-»TtlH  .U'11  • 

n,  CONTINUE 

p%  TS‘.-P-f' 

yHITE  (f  •“>  , p.  , , u .1.  1 

N f 0P>“  I ' ' ‘ , 

,btte 

p f OPN*  T I / l»  • 


24 
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l* 

sus»nuT  tnf  7 7 is 

2‘ 

c 

Xm 

c 

This  P»nGR»"  Cfirt 

«• 

c 

s* 

INTFGFR  ».u 

OIMFNSION  * 1*  7S 

?• 

ht : (s  - 1 1 • r«- 1 ) 

• • 

c 

c 

GENr»4Tr  L -f'UNC  T TON  ' 

10* 

c 

ll« 

DO  700  jrl.S 

12- 

* (J. J1 ■ 1 

n* 

500 

COSTISUF 

J - s*  1 

15» 

«TT:^.«T-I 

IK* 

NK  : N.| 

I»* 

c 

!•* 

c 

smf»»T'-  L-ruNCTIOH  1 

H* 

Ti*~ 

c 

• 0 *01  H >rt  .Nfe 

21* 

DO  10  H3j,»TT 

22- 

« (M.M  1 ) = 1 

2^. 

n 

f ONT  TNUF 

2»* 

K2=Kt«l 

2S* 

►4-J 

26* 

DO  1 1 T:k  ? .N 

2 7* 

z *«  . 1 

21* 

no  1?  1 1-- 1 

2^. 

A |H  , I 1 r T I 

I»* 

c 

M* 

c 

COUNT  NUNUr*  or  SF  T<; 

52* 

c 

«I* 

h:m*  1 

»•* 

12 

CONT  JNUF 

IS* 

1 I 

rONT  TNUF 

SK* 

Jt  M 

! 7* 

MTT=H-*1  ♦mF«|N>i 

SI* 

STI 

CONT  IN'JF 

7S* 

H Z M.l 

*t* 

»f TU»N 

•1* 

CND 
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MOD-ANOVA:  An  Analysis  of  Variance  Program 

for  2*^  and  3^  Factorial  Experiments 
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